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The spatial structure of the axonal and dendritic arborizations is closely related to the func- 
tionality of specific neurons or neuronal subsystems. The present work describes how multiscale 
Minkowski functionals can be used in order to characterize and compare the spatial organization 
of these two types of arborizations. The discrimination potential of the method is illustrated with 
respect to three classes of cortical neurons. 
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Human intelligence is to a great extent determined by 
the connectivity between the myriad of neuronal cells in 
the nervous system (Cajal:1989). Biologically, these con- 
nections are implemented by the synapses between neu- 
rons, which are strongly related to the shape of axons 
and dendrites (e.g. [UEIE]). For instance, neural cells 
exhibiting more intricate dendritic or axonal arboriza- 
tions tend to promote more connections. Therefore, the 
systematic study of the morphological properties of the 
axonal and dendritic trees can provide essential informa- 
tion about the connectivity of the nervous system. Be- 
cause neurite outgrowth can be understood as a complex 
dynamical system involving pattern formation, biochem- 
ical regulation and electrical activity, they are particu- 
larly suitable for characterization and modeling by using 
concepts and methods from physics, computation and bi- 
ology. 

As the spatial distribution of synapses is to a great ex- 
tent influenced by the interaction between the shape of 
dendrites and axons, and also because they share several 
developmental aspects, one question of particular rele- 
vance concerns the study to which extent these two types 
of structures arc related (i.e. similar or distinct). Inter- 
estingly, both axons and dendrites start their lives as neu- 
rites, being differentiated during development. Along the 
early stages of neuronal growth, which involve many out- 
growths and retractions, one of these neurites dominates 
and become the cell axon, assuming distinct biomolec- 
ular composition and function. In particular, axons 
seek actively for targets while being influenced by many 
trophic factors including differential adhesiveness, gal- 
vanotropism and chemotropism (e.g. |H El E]). Conse- 
quently, the shape of axons is directly affected by external 
field influences. At the same time, both dendrites and ax- 
ons share intrinsic genetic and biochemical basis, which 
constrains neuronal growth and may imply intrinsic sim- 
ilarities between the dendritic and axonal arborizations 
of a same cell or category of cells. The current arti- 
cle addresses this important problem from the objective 
perspective of neuronal shape quantification by additive 
functionals [7\, which have been previously applied with 
encouraging success to the characterization of dendritic 



morphology (e.g. [5JE]). 

Integral geometry algorithms have been successfully 
used to characterize morphologically complex patterns 
and structures whose process of formation is not precisely 
known and is subject of modeling [TJ. The central proce- 
dure is the calculation of intrinsic volumes or Minkowski 
functionals, a generalization of the usual determination 
of volume. 

The Minkowski functionals of a body K in the plane 
are proportional to familiar geometric quantities, namely 
its area A(K), perimeter U(K) and the connectivity or 
Eulcr number x{K)- The usual definition of the connec- 
tivity from algebraic topology in two dimensions is the 
difference between the number of connected n c compo- 
nents and the number of holes n^, x(-^0 = n c~ n h- I n the 
Euclidean space, there is an additional geometric quan- 
tity, the mean curvature or breadth. Moreover there are 
two kinds of holes to consider: A pure hole, which is a 
completely closed region of white voxels surrounded by 
black pixels and tunnels. The Euler characteristic for the 
Euclidean space is then given by 

X{K) = n c - n t + n h , (1) 

where n t is the number of tunnels and nh is the number 
of pure holes. One instance where these functionals ap- 
pear naturally is while attempting to describe the change 
in volume as the body K, now assumed to be convex, un- 
dergoes a dilation through a parallel set process using a 
ball B r of radius r 

47T 

V(K 8 B r ) = V(K) + S{K)r + 2nB(K)r 2 + — r 3 . (2) 

o 

Generalizing to higher dimensions, the change in hyper 
volume is given by the Steiner formula 

v d (K®B r ) = Y J ( d )wl d \K)r», (3) 
where the coefficients Wu are referred to as Minkowski 
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functionals. For instance in the Euclidean space (d = 3), 

W^(K) = V(K), wi 3) (K) = S ^l, (4) 
W^(K) = ^B(K), W^(K) = ^- X (K). 

Despite the wealth of results and continuum formulae 
for obtaining these functionals, it is useful to resource 
to the discrete nature of the binary images we wish to 
analyze by looking at the distribution of voxels in a Eu- 
clidean spatial lattice. By exploring the additivity of 
the Minkowski functionals their estimation resumes to 
counting the multiplicity of basic building blocks that 
disjointedly compose the object. The fundamental infor- 
mation needed here is a relationship for the functionals 
of an open interior of a n-dimensional body K which is 
embedded into a d-dimensional space 

W { u d] {k) = (-l) d+n+v W^(K),v = 0, . . . , d. (5) 

With the absence of overlap between these building 
blocks and using the property of additivity of these func- 
tionals, we may write for a pattern V composed of dis- 
joint convex interior pieces N m , 

Wi J d \V) = Y,W^ d) {N m )n m {V), u = 0,...,d , (6) 

I a 

where n m {T') stands for the number of building elements 
of each type m occurring in the pattern V. For the 
three-dimensional space we display in Table |T] the value 
of Minkowski functionals for the building elements in a 
orthogonal lattice of voxels. By using the information 
(with a = 1) presented in Table [I] and Equations Q, ^ 
and (JsJ we have 

V = n 3 , S=-6n 3 + 2n 2 , (7) 
2B = 3n 3 - 2n 2 + Hi, X = ~ n 3 + n 2 -ni+ n Q , 

Where n 3 is the number of interior cubes, n 2 is the num- 
ber of open faces, n\ is the number of sides and no is 
the number of vertices. So the procedure to calculate 
Minkowski functionals of a pattern V has been reduced 
to the proper counting of the number of elementary bod- 
ies of each type that compose a voxel (cubes, faces, edges 
and vertexes) involved in the make up of V . 
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TABLE I: Minkowski functionals Mi' for open elementary 
open sets N m which compose a voxel K(of side a — 1): V 
(vertex), L (side), F (face) e C (cube). 



The recent effort aimed at replicating a working micro 
column of cortex tissue |10j has motivated the dissemi- 
nation of comprehensive data bases containing different 
anatomical classes of neuronal cells. We have consid- 
ered several neurolucida data from anatomically differ- 
ent types of cells, namely the anatomical classes MCI, 
BTC and LBC that populate the cortex and are engaged 
in distinct functions in the working of a neocortical mi- 
cro column. The morphological data was obtained from 
the Brain Mind Institute, EPFL database [IT]. Exact 
dilations of both the axonal and dendritic trees were 
performed until the majority of holes and tunnels dis- 
appeared, i.e., when the Euler number approaches 1. 
This implied a total number of 1000 dilations. All four 
Minkowski functionals were calculated along such dila- 
tions, so as to obtain a morphological signature for the 
dendritic and axonal arborizations of each cell. We then 
conducted a morphological analysis using a subset those 
values at selected scales, namely after 1,250, 500, 750 and 
1000 exact dilation steps. In order to reduce the effect of 
the size of the cells on the measurements, the function- 
als were normalized respectively to the adequate power 
of the diameter of the cell, i.e the volume functional was 
normalized by I? , the surface area was normalized by L 2 
and the mean-curvature was normalized by L. The Euler- 
number is a dimcnsionlcss topological measure and con- 
sequently did not undergo normalization. All measures 
were subsequently standardized, which is accomplished 
by subtracting the average and dividing by the standard 
deviation of each respective type of measurement (i.e. 
one of the 4 functionals). The standardized values there- 
fore have zero means and unit variance, and most of their 
values are comprised between -2 and 2. Such a stan- 
dardization reduces the influence of the overall relative 
magnitude of the different measurements. 

The functionals were calculated separately for the den- 
dritic and axonal trees of each neuron. For each case, the 
obtained measurements were organized into morphologi- 
cal vectors in the 20-dimensional space defined by the five 
spatial scales of each of the four functionals. Therefore, 
given a neuronal cell represented by their respective mea- 
surement vector, the morphological difference between 
its axon and dendritic arborizations can be quantified in 
terms of the magnitude [i of the vector corresponding 
to the difference between the respective feature vectors 
(i.e. axon and dendrites). Another important descriptor 
of the morphological space which has been considered 
in this work is the angle (f> between the above difference 
vector and the vector which represents the mean of all 
such difference vectors. Such a measurement provides 
and indication of how much the axon-dendrite in each 
cell departs from the overall prototypical case (i.e. the 
average) . 

The result of the above procedure is presented in the 
scatterplot (i.e. the angle versus the magnitude of differ- 
ence vectors for each neuron) in the composited Figure [T| 
which also shows the relative densities of the magnitudes 
of the difference vectors obtained by using a simple Gaus- 
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sian fit. Recall that total morphologic similarity between 
the axonal and dendritic trees for each cell would result 
in a magnitude of difference vector equal to zero. The 
scatterplot in Figure [T| shows that most cells resulted 
with magnitude between and 2 and angle between 
and 6, with an outlier with angle near to 10 in the case 
of the MC2 cell category. As shown by the respective 
densities, the three classes of neuronal cells resulted with 
distinct dispersions, with the MC2 class being character- 
ized by the smallest variance while the BTC class shows 
the largest variance. In order to obtain a more objective 
indication of the separation between the three neuronal 
cell categories, we performed the ANOVA variance test 
(e.g. [T5]) and obtained the results (p— values) shown in 
Table Hfl These values indicate that the MC2 and LBC 
classes are those less likely to have come from the same 
population of cells, with a respective p— value of only 
0.0088. Substantially higher values were obtained for the 
other pairs of categories, suggesting less separability in 
those cases. 

Figure [T] also shows an example of neuronal cell for 
each of the anatomical classes, lying close to the center 
of the respective Gaussian distribution fit. Despite the 
difficulty of comparing the original 3D axonal and den- 
dritic arborizations from such 2D projections, the larger 
difference between the two arborizations obtained in the 
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TABLE II: The resulting p- values obtained through ANOVA 
for the 3 classes of neurons. 



case of the LBC category is still clear, as well as the 
greater similarity observed for the MC2 cell example. 

The obtained results suggest varying degrees of mor- 
phological relationship between the axonal and dendritic 
arborizations among all the considered individual cells 
and also between the three classes BTC, LBC and MC2. 
In the latter case, the classes LBC e MC2 resulted as be- 
ing strongly statistically distinct. In brief, the method- 
ology reported in this work has been found to provide 
a sensitive means for inferring differences between the 
axonal and dendritic structures. This is a particularly 
important result, as such differences may be related 
to important differentiation factors during the neuronal 
growth, including those caused by genetic modifications 
(cell predestination), environmental constraints (e.g., in- 
teractions with other cells and target specificity) , as well 



4 



as the history of presented stimuli. The identification 
of such differences also present potential comple- 
mentary subsidy for neuronal cell classification and di- 
agnosis of abnormalities in neurological diseases, which 
correspond to interesting topics for further research. 
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